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Abstract 

We examine the dependence of parton distribution functions (PDFs) on the value of the QCD 
coupling strength a s (Mz). We explain a simple method that is rigorously valid in the quadratic 
approximation normally applied in PDF fitting, and fully reproduces the correlated dependence of 
theoretical cross sections on a s and PDF parameters. This method is based on a statistical relation 
that allows one to add the uncertainty produced by a s , computed with some special PDF sets, 
in quadrature with the PDF uncertainty obtained for the fixed a s value (such as the CTEQ6.6 
PDF set). A series of four CTEQ6.6AS PDFs realizing this approach, for a s values in the interval 
0.116 < a s (Mz) < 0.120, is presented. Using these PDFs, the combined a s and PDF uncertainty 
is assessed for theoretical predictions at the Fermilab Tevatron and Large Hadron Collider. 
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I. INTRODUCTION 



The global analysis of quantum chromodynamics (QCD) refers to the use of data from 
many short-distance scattering processes to construct, within some approximations, universal 
parton distribution functions (PDFs) that can be used to calculate hadronic cross sections 
for the QCD and electroweak theories. An important part of the analysis is to determine the 
uncertainties of the PDFs. The uncertainties have several sources, both experimental and 
theoretical. An example of an uncorrelated experimental uncertainty, i.e., an uncertainty 
induced by data, is the statistical error on the cross section measurements used as the input 
to the global analysis. An example of a correlated experimental uncertainty, also induced by 
data, is the luminosity error, which modifies the normalization of all cross sections measured 
by one experiment by the same factor. An example of a theoretical uncertainty is the choice 
of a momentum scale Q used in the theoretical calculations for some process; this uncertainty 
depends on the order of perturbative approximation (next-to-leading order, NLO, in a s for 
the CTEQ analysis in this paper) and is highly correlated between different points in the 
phase space of the process. Methods for calculating these various uncertainties have been 
developed for the CTEQ analysis [if] and other analyses of PDFs 0-0] • 

Two overall theoretical uncertainties of the global analysis are the "parametrization error" 
and the uncertainty in the QCD coupling strength a s . The first concerns the parametric form 
for the nonperturbative PDFs at a low momentum scale Qq, used as a boundary condition 
for predicting the PDFs at other scales through their perturbative evolution. The choice of 
the functional form affects the predictions and so contributes an uncertainty to the published 
PDFs. 

The strong coupling a s is another basic parameter that affects every theoretical cross 
section in the global fit. As with the PDF parametrizations, one must decide how to evaluate 
the scale dependence of the running coupling constant a s (Q) and provide a value of a s at 
some scale Q' as a boundary condition. Both choices contribute additional uncertainties 
to the determination of PDFs. The Q dependence of a s (Q) is not uniquely determined, 
but estimated from the NLO approximation of the renormalization group (RG) equation for 
a s (Q): 

da s O 2 /?i 3 , 4 

dQ AH 07T Z 

where (3 = 11 — (2/3)n/, and (3\ = 102 — (38/3)n/. In NLO perturbation theory, we neglect 
terms of order af. Then there are different solutions of ([I]), formally equivalent in the NLO 
approximation, but differing in higher orders of perturbation theory. Previous studies have 
shown that the associated uncertainty, coming from the approximation for the Q dependence 
of a s (Q), is small compared to other sources of the PDF uncertainty [?|]. 

In this paper we are concerned with the value of a s at a typical scale (Q) for hard 
interactions. The differential equation ([1]) requires an initial value, commonly taken from a 
combination of precise experimental measurements. The conventional choice is the value of 
a s (Q) at the Z boson mass, Q = Mz- Any value of Q could be taken to set the initial value, 
but Q = Mz is the natural choice, because Mz is known precisely, and some experiments 
determine the strong coupling directly at this momentum scale. 

The input value of a s (Mz) and its uncertainty are continually being refined, as new 
experimental data are obtained. The global PDF analysis itself could "determine" (i.e., 
estimate within the errors) the value of a s (Mz) from the comprehensive hadronic data it 
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examines. However, the resulting determination turns out to be quite uncertain: the included 
scattering processes (deep-inelastic scattering, Drell-Yan process, pp — > jets, etc.) have 
limited accuracy for measuring a s (Mz), precisely because they depend on uncertain PDFs. 

The most accurate value of a s (Mz) is instead available from a worldwide compilation 
provided by the Particle Data Group (PDG) jsl, The 2009 world average value of a s 
based on eight different measurement techniques is [8[ 

a s (M z ) = 0.1184 ± 0.0007. (2) 

Most experimental data that are an input to this world average are independent of parton 
distributions. For example, precision data from e + e~ annihilation at LEP provide strong 
and direct constraints on a s (Mz), independently of the PDFs. The uncertainty range in 
Eq. 02]) nominally corresponds to the 68% confidence level (CL); yet, when assigning this 
level, one must be aware that the analysis mixes NLO and NNLO calculations in various 
techniques and relies on some theoretical models. Even with these caveats, the accuracy of 
the world-average value is clearly superior to the value determined solely from the hadronic 
scattering data. 

Traditionally, our analysis presents the best-fit PDFs and their parametrization uncertain- 
ties for a constant value of a s (Mz) close to its latest world-average central value. Separately, 
the uncertainty in the PDFs induced by the uncertainty in a s {Mz) is assessed, by producing 
a few alternative PDF fits for a range of values of a s (M z ) around the central a s (M z ) value 
0]. We add the two uncertainties in quadrature to estimate the total uncertainty result- 
ing both from the PDF parametrization and a s . In this paper, we examine this procedure 
(designated as "fitting method 1") in detail, in order to determine how well it captures the 
correlation, or interplay, of the a s uncertainty and PDF uncertainty. 

On the other hand, the most complete method to evaluate the combined PDF + a s 
uncertainty in the global fit is to vary the theoretical value of a s [Mz) as an additional 
fitting parameter, while including the world-average value of a s (Mz) with its experimental 
uncertainty as a precise experimental constraint on a s (Mz), in combination with the rest 
of the hadronic data. The uncertainty from a s can then be deduced by this technique - 
designated as "fitting method 2" - on the same footing, and by the same Hessian techniques, 
as the uncertainties in the PDF parameters. [l| 

Our main finding is that the simpler framework of the fitting method 1 (based on the ct s 
series) is sufficient to reproduce all correlated dependence of a s and PDF parameters arising 
in the fuller treatment of the fitting method 2. Consequently, the a s series of CTEQ6.6 fits, 
generated in fitting method 1, is sufficient for determining the total PDF+a s uncertainty in 
theory calculations, e.g., predictions for the Large Hadron Collider (LHC). This conclusion 
is borne out both numerically and formally, within the quadratic approximation for the log- 
likelihood function \ 2 i n the vicinity of the best fit. The total PDF + a s uncertaint y ca n be 
thus calculated without resorting to the alternative, quite elaborate, methods (1, [Hi U|, but 
simply by computing the theoretical cross sections for two additional CTEQ6.6 PDF sets for 
the extreme values of a s , besides the usual cross sections for the 44 CTEQ6.6 eigenvector 
PDFs. 

Our conclusion that the addition in quadrature is sufficient, regardless of the magnitude 
of the correlation, follows from an easily overlooked distinction between the "absence of 
correlation", on one hand, and "the independent addition of the a s uncertainty and PDF 
uncertainty" , on the other hand. As we show, variations in ct s generally induce compensating 
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adjustments in the preferred PDF parameters (correlation) to preserve agreement with those 
experimental data sets that simultaneously constrain a s and the PDFs. At the same time, it 
is possible to define an u a s uncertainty" that quantifies all correlation effects, is independent 
from the PDF uncertainty for the best-fit a s , and can be added in quadrature. 

Section [III presents the new a s series of PDFs, based on the 2009 world average value 
of a s (M z ) = 0.118 and four additional PDF sets found for 0.116 < a s < 0.120. This 
series, denoted as CTEQ6.6AS, extends the CTEQ6.6 series of PDFs lUl; it uses the same 
input data and methods as the CTEQ6.6 study. In the same section, we also introduce an 
alternative fit, designated CTEQ6.6FAS, with a floating a s , constrained by the world-average 
value as an experimental input. We provide a heuristic explanation of why the CTEQ6.6AS 
and CTEQ6.6FAS methods should result in close estimates of the total PDF + a s uncertainty. 
The detailed proof is given in the Appendix. In Section IHfl we compare the CTEQ6.6AS 
and CTEQ6.6FAS methods numerically. Also, applications based on the resulting PDFs are 
shown, including our estimates for the current a s uncertainty on select theoretical predictions 
for the LHC. 

Section IIVI addresses a separate issue, the degree of the correlation between the a s and 
PDF uncertainties imposed by the hadronic data in the current fits. We elucidate this issue 
by applying the correlation analysis of Refs. [l|, Qjl, 13] to explore which data sets in the 



CTEQ6.6 sample impose the tightest constraints on a s . The paper concludes in Section IVl 
with our recommendation to use the a s series, CTEQ6.6AS, combined in quadrature with 
CTEQ6.6, to assess the effect of a s uncertainty on theoretical predictions. 



II. QCD COUPLING AS AN INPUT AND OUTPUT OF THE PDF ANALYSIS 
A. Diagonalization of the error matrix with respect to a s 

Perturbative QCD cross sections calculated to next-to-leading order (NLO) have un- 
certainties arising from (1) the residual scale dependence of the cross section due to the 
truncation of the series at NLO, (2) the PDF uncertainties derived from the global PDF 
fits and (3) the uncertainty in the value of a s (Mz) used in the cross section calculation. 
The residual scale dependence will depend on the individual perturbative cross section being 
evaluated. The PDF uncertainty can be determined using the PDF eigenvector sets provided 
by the global fitting groups. CTEQ and MSTW provide eigenvector PDFs (44 for CTEQ6.6 



12J and 40 for MSTW2008 [10J) determined using the Hessian method, while the NNPDF 
collaboration provides ~ 1,000 PDFs|§|. The value of a s (Mz) used in any cross section 
evaluation must be the same as its value in the global fit providing the PDFs. Thus, the 
uncertainty in a s (Mz) must be evaluated by using the PDF sets in which the same value of 
a s (M z ) has been assumed. 

The treatment of the a s uncertainties in any NLO cross section evaluation depends on 
their relative size compared to the intrinsic PDF uncertainties, and on whether there exists 
a correlation between the value of a s [Mz) and the Hessian error PDFs. One cannot exclude 
the possibility that the PDF and a s uncertainties affect each other. For example, it has been 
known for a long time that the shape of the gluon PDF found from the global fit changes 
considerably if a s (Mz) is changed in the fit. 

There are two prevailing approaches to the choice of the value of a s (Mz) in the global 
PDF fits: either the world average is taken as an input, or a s (Mz) is determined as an output 
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of the global fit. The CTEQ, HERAPDF and NNPDF sets of PDFs use the first approach, 
while MSTW uses the second. Related to this choice is the evaluation of the uncertainty 
in a s [Mz)'- either the uncertainty can be directly related to that determined by the world 
average, or the uncertainty can be determined by examining the impact of varying a s [Mz) 
in the global fit. Again, CTEQ and the NNPDF collaborations have used the first approach, 
and MSTW the second. 

The reason behind the first approach is that the most precise processes contributing to 
the world average (r and quarkonium decays and e + e~ event shapes at the Z-pole energy) 
are free of the PDF uncertainties and place stronger constraints on the value of a s (Mz) than 
the hadronic scattering data alone. Knowledge gained from those processes, summarized 
in the most recent PDG value of a s (Mz), naturally provides a useful external input to the 
global QCD analysis. After all, the idea of global analysis is that QCD is a fundamental 
theory that describes all aspects of strong interactions. The value of a s {Mz) is universal. 
The value that we use should be the most accurately determined value. 

The second approach has the advantage of being independent from other sources. It 
circumvents the issue of combining constraints on a s (Mz) from heterogeneous measurements 
and mixing different theoretical frameworks in the combination procedure. The price to pay, 
however, is the loss of the constraining power supplied by the most precise measurements. 

As an example of the uncertainty in a s (Mz), consider again the 2009 world average [§[ 
in Eq. j2]), its determination is dominated by the processes that do not require the PDFs. 
The interpretation of the confidence level assigned to this uncertainty (nominally 68%) is 
still not fully settled for the reasons mentioned in Section IB For the purposes of this study, 
we adopt a somewhat more conservative estimate of the uncertainty proposed at the 2009 
Les Houches workshop [lH and adopted by the PDF4LHC working group flo| . 

a s (M z ) = 0.118 ± 0.002 (90% CL), (3) 

corresponding to an uncertainty of ±0.0012 at the 68% CL. 

Compare that to the value obtained by MSTW from fitting a s in the PDF analysis 0, 10 



a s (M z )= 0.1202 i +°;°^ (68% CL), (4) 

a s {M z ) = 0.1202 j (90% CL) (5) 

at next-to-leading order; or 

a s (M z ) = 0.1171 ± 0.0014 (68% CL) (6) 

at next-to-next-to- leading order. Or, compare that to our own determination of a s (Mz) = 
0.118 ± 0.005 (90% CL) from a fit without the world-average constraint that is described 
below, cf. Eq. ([TBI . The result is consistent with that in Eq. ([2]); but the uncertainty on 
a s (Mz) is considerably larger if a s is fitted without the world-average constraint. The central 
value of a s [Mz) returned by this fit is 0.118, which coincides with the fixed value used in 
the CTEQ6.6 PDF analysis. 

One can also envision a third, most general approach, in which the a s (Mz) range measured 
by the most precise techniques, i.e. the world average value, is included as an input together 
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with the usual hadronic scattering data; and the theoretical value of a s (Mz), fitted as a free 
parameter together with the PDF parameters, is returned as an output constrained by the 
combination of all measurements. Such a fit is the most direct in quantifying improvements in 
the accuracy of the output a s [Mz) resulting from the hadronic scattering data (as compared 
to the more precise constraints); as well as in probing the correlation between a s [Mz) and 
PDF parameters. 

In this paper, we examine the CTEQ global hadronic data and the world-average a s (Mz) 
in Eq. according to this more general approach, and apply the usual Hessian technique |l| 
to study the combined PDF + a s uncertainty (i.e. the usual PDF shape parameters plus 1 
more parameter for a s ). However, when implemented straightforwardly, this approach runs 
into a practical inconvenience: each of the eigenvector PDFs is associated with its own value 
of a s {M z ). 

The quadratic approximation provides a remarkable bypass for this shortcoming. By 
applying the Data Set Diagonalization method introduced in Ref. [13], it is possible to re- 
diagonalize the parameter space so that only one eigenvector corresponds to the change of a s , 



whereas all the other eigenvectors are immune to the change in a s . (As shown in Ref. 17J , this 
conclusion also holds for any other distinct parameter in the global analysis.) In other words, 
it is easy to construct a pair of additional PDF sets corresponding to the maximal tolerated 
excursions of a s . The a s uncertainty computed from the difference of these PDFs can be 
added in quadrature to the CTEQ6.6 uncertainty to reproduce the combined uncertainty 
with the full PDF-a s correlation. We refer to the reduced computation as "fitting method 
1," and the full computation as "fitting method 2." 

A formal proof that methods 1 and 2 are equivalent, within the quadratic approximation, 
is given in the Appendix. Here we provide a heuristic argument that illustrates the choice 
of the two PDF eigenvector sets probing the a s uncertainty. 

The a s parameter in the full fit is associated with a new twenty-third direction that is 
orthogonal to the hyperplane of the PDF parameters for the central value of a s {Mz) = 0.118; 
i.e., the 22-parameter space probed by the CTEQ6.6 PDF set. Originally, the orthonormal 
eigenvector PDFs of the CTEQ6.6 fit correspond to a representation such that the excursion 
by the same distance from the best fit in the 22-dimensional hyperspace results in the same 
increase of x 2 \ that is, the surfaces of constant x 2 are 22-dimensional hyperspheres. 

If a s is allowed to deviate from its best-fit value, a preferred direction emerges in the 
22-dimensional hyperspace, along which A% 2 grows most slowly, because changes in a s are 
compensated by changes in some linear combination of the PDF parameters. Rotate the 
CTEQ6.6 basis so that the "physical" PDF parameter a\ corresponds to the variation of 
this preferred combination from its best-fit value, and analogously, the "physical" parameter 
cio corresponds to the variation of ct s (Mz) from its best-fit value (which remains at 0.118). 
With this choice, ao = and a± = at the best fit. All correlations between ct s and 
PDFs are encapsulated in the x 2 dependence on ao and a±. Let us focus on this dependence. 
Uncertainties due to the other (combinations of) PDF parameters 02, 022 are independent 
of a s and can be added in quadrature at the end. 

In the {a , ai} plane, the set of the allowed PDFs corresponds to the inside of an ellipse 
< T 2 , where T is the tolerance parameter or some other parameter defining the typical 
X 2 at the boundary of the allowed region. The CTEQ6.6 fit is equivalent to probing the 
inside of the ellipse for a fixed a Q = 0, as shown in the left inset of Fig.[TJ The contribution 
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Physical basis a* 

Ax 2 = J2i,j Hijdidj 



ai d\ 




AI? = \ (X(B) - X(D)f AI 2 = \ (X(A) - X(C)f 

FIG. 1: Imagine fitting two free parameters. The left figure shows uncertainties computed in a fit 
with oq = and free a\. The right figure shows a x 2 scan over ao with free a\. 



Physical basis a; 

Ax 2 = Eij H i}j didj 
ai 



Eigenvector bases y it z\ 





Ax 2 
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- X(C)) 2 + - X(D)) 2 

= AXg + AlJ 



FIG. 2: Imagine fitting two free parameters. Relations are illustrated, between the basis of physical 
parameters {ao,ai}, the principal axis basis {yo,yi}, and the a s -excursion basis {zq,z\}. 
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of the ai direction to the CTEQ6.6 uncertainty for an observable X can be evaluated as 

AXl = \{X{B)-X{D))\ (7) 

where X(B) and X(D) are the values of X at points B and D in the figure. 

Next, consider the minimal and maximal excursions of the parameter ao allowed inside 
the ellipse, for a free parameter a%. These excursions are reached at points A and C shown 
in the right inset of Fig.HJ The uncertainty along this direction is 

AX 2 Q = 1 -(X(A)-X(C)f. (8) 

In practice, points A and C are found by a scan of the dependence of Ax 2 (ao, a\) on ao, 
for a varying a±. To describe the uncertainty of an observable X, we rotate the {a , a±} basis 
to a basis {yo, y{\ of eigenvectors of the Hessian matrix. The orthogonal directions yo and yi 
specify the principal axes of the ellipse, as seen in the left inset of Fig.[2J Since the maximal 
range of the ao parameter is specified by the line segment AC in Fig. Q] or Fig-El we could 
further rotate the {2/0,2/1} basis to another basis {,20,-21}, after rescaling the eigenvectors yo 
and yi to unity, as shown in the right inset of Fig.[2j The vector zq is chosen to be along the 
AC direction. It is shown in the Appendix that the vector z\ is along the direction of the 
line segment BD. It is perpendicular to the zq direction in the {zo,z\} representation, even 
though the z and zi directions are not perpendicular in the {a , a{\ representation. Taking 
this conclusion for granted here, we could easily see from the left inset of Fig.[2]that there is 
no a dependence along the z\ direction. From the right inset, the total uncertainty of the 
observable X can be computed as 

AX 2 = \ [(X(A) - X{C)f + (X(B) - X{D)f] (9) 

= AX^ + AXl (10) 

This uncertainty is equal to the quadrature sum of the uncertainties along the AC and BD 
directions, as has been stated. 

We see that the zo and zi directions, which are initially not orthogonal in the {ao,ai} 
basis, can be made such by rotation and scaling. Alternatively, the orthonormal {zq, z{\ basis 
can be obtained from the {ao, ai} basis by a shear transformation in the negative a\ direction, 
followed by a scaling transformation along the ao direction. In the 23-dimensional case, the 
condition a = z = defines the hyperplane spanned by the 22 CTEQ6.6 PDF parameters, 
for the fixed best-fit a s . This hyperplane is made orthogonal to z by a rotation and scaling 
transformation, or, equivalently, by a shear and scaling transformation, analogously to the 
2-dimensional case. See the Appendix for the full discussion. 

B. Fits with a variable a s : explicit realizations 

We will now explicitly construct two next-to-leading order (NLO) fits of the kinds de- 
scribed above, using the CTEQ6.6 sample (l2| of hadronic data 1 . To this end, we modify 



1 Like in all CTEQ global fits, we include full NLO matrix elements in DIS and most Drell-Yan observables. 
In calculations for inclusive jet production and W lepton asymmetry, where the full NLO/resummed results 
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the setup of the CTEQ6.6 analysis to allow a s (M z ) to vary within the global fit, and to con- 
strain these variations by the world- average (w.a.) value (oc s ) ± (5a s ) w = 0.118 ± 0.002 
(at 90% CL) included as a separate data input, in addition to the complete CTEQ6.6 set of 
hadronic scattering data. Agreement with this precision data value is just as desirable (or 
more so) as agreement with individual data points in the hadronic data sets. So, to assure 
this agreement, we add a new contribution Xa to the global log-likelihood function in the 
global analysis: 

X 2 = X 2 CTEQ6.6 + ^' ( n ) 

where 

a s (M z ) - MM z )) wa ] 2 

The Xa s t erm is multiplied by a weighting factor A to match the confidence interval of 
the world-average a s with the tolerance on the increase in x 2 allowed for acceptable fits. 
We choose A so that a s (M z ) values outside of the 90% CL interval, a s (Mz) < 0.116 or 
a s (Mz) > 0.120, result in a penalty beyond the tolerance for the increase in x 2 ■ 

Without the contribution of the world-average value, the constraints of the global fit on 
a s are relatively weak, 

a s (M z ) = 0.118 ± 0.005 (90% CL). (13) 

The central value returned by this fit is practically identical to either the PDG or Les Houches 
workshop central values. When the world-average constraint is included, the final result of 
the global analysis changes to 

a s (M z ) = 0.1180 ± 0.0019 (90% CL). (14) 

Again, its central value is practically the same, but the uncertainty is much smaller. Thus 
the constraint on a s (Mz) is dominated by the world-average uncertainty, (Sa s ) w _ a _. 

We explore the vicinity of the best fit in two ways. First, we construct best-fit PDFs sets 
for four alternative values of a s (Mz), 

a s {M z ) = 0.116, 0.117, 0.119, and 0.120. (15) 

These PDFs are named as 

CTEQ6.6AS = AS„ 2l AS- h AS +1 , and AS +2 . (16) 

The CTEQ6.6AS PDFs for the two extreme variations, AS-2 and AS+2, correspond 
to slightly more than two standard deviations according to the PDG error in Eq. fl2]), or 

are prohibitively CPU-extensive, K-factor tables are used to look up the ratio of the NLO cross section 
to the LO cross section (with the LO cross section calculated using the NLO PDFs, and using the 2-loop 
a s ) separately for each data point. The look-up tables depend very weakly on the input PDF parameters. 
They are updated in the course of the fitting to ensure that the K-factors have not drifted from their 
initial values, and we again check that the calculations retain their full NLO accuracy at the end of the 
global fit. 



xl 3 = A 
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approximately to the 90% CL uncertainty according to the 2009 Les Houches prescription 
in Eq. Q. The intermediate PDF sets, AS_i and AS + i, provide additional information on 
the a s dependence. 

Alternatively, the combined PDF and a s uncertainty is quantified by the diagonalization 
of the Hessian matrix, in terms of 46 extreme PDF eigenvector sets for 23 independent 
combinations of theoretical parameters. Each eigenvector set is associated with its own 
a s (Mz) value within the best-fit range of Eq. (Tl4|) . This series of eigenvector PDFs is called 
CTEQ6.6FAS. It will be compared with the CTEQ6.6AS series in the next section. 

III. NUMERICAL RESULTS 

A. Comparison of the PDF and a s uncertainties 

The AS_2 and AS + 2 PDFs of the CTEQ6.6AS series for the gluon, u quark, and s quark 
are compared with the CTEQ6.6 PDF uncertainty band for the momentum scale Q = 2 
GeV in Fig. [31 and for Q = 85 GeV in Fig.|H Each figure shows the ratio of f(x, Q) in the 
AS n set to the CTEQ6.6M PDF of the same flavor (solid and dashed curves), as well as the 
asymmetric CTEQ6.6 PDF uncertainty for this flavor (shaded region), as functions of the 
momentum fraction x. The upper and lower boundaries 1 ± [Af(x, Q)]±/f{x, Q) of the PDF 
uncertainty band are given by the asymmetric PDF errors, as 



22 

]T[/^,Q)-/o(x,Q)] 2 for fi>f , 
i=i 

22 

E l/fo Q) ~ /<>(*> Q)f for fi < ^ (17) 

i=l 

For the shown partons (g, u, and s), the PDF uncertainty is larger than the a s dependence 
over the range 0.116 < a s (Mz) < 0.120. The figures for the d quark (not shown) looks similar 
to the figures for the u quark. The relation between the PDF and a s uncertainties for u 
and d quarks is in between these relations for the u and s quarks. Finally, the c and b 
uncertainties are qualitatively similar to those for g. 

We now add the CTEQ6.6 PDF uncertainty in quadrature to the CTEQ6.6AS a s uncer- 
tainty, and compare this sum to the uncertainty of the 23-parameter fit CTEQ6.6FAS with 
the floating a s and the world average constraint. The result is shown in Fig.[5j The figure 
compares the PDF uncertainty bands from the two realizations of the a s fit, for the parton 
flavors that show the largest a s dependence: gluons and charm quarks. The CTEQ6.6 er- 
ror band is indicated by dashed lines; while the two realizations of the a s fit by a hatched 
error band with solid borders and a filled error band with dotted borders. We observe only 
small differences between the CTEQ6.6+CTEQ6.6AS and CTEQ6.6FAS PDFs for the fla- 
vors shown in the figure. The differences are even smaller for the other quark PDFs. Thus, 
the two realizations of the a s series produce nearly identical results, confirming the adequacy 
of adding the PDF and a s uncertainties in quadrature. 



[A/(x,Q)]. 
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[A/(z,Q)L 
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g at Q=2 GeV 




FIG. 3: Comparison of the PDF and a s uncertainties for the gluon, u quark, and s quark, at 
Q = 2 GeV. The vertical axis is the ratio of f(x,Q) to the CTEQ6.6M best fit. The CTEQ6.6 
PDF uncertainty range for the central a s value, 0.118, is shown as the shaded region. The ratio of 
AS+2 and AS-2 PDFs to the corresponding CTEQ6.6M PDFs are shown as the solid and dashed 
curves, respectively. 
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FIG. 5: Comparison of CTEQ6.6+CTEQ6.6AS (blue solid) and CTEQ6.6FAS (green hatched) 
uncertainty bands and CTEQ6.6 uncertainty bands (indicated by black dashed lines), normalized 
to the standard CTEQ6.6M fit, for gluon and charm PDFs at Q = 2 and 85 GeV. 
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B. Uncertainties of cross section predictions 



Using the a s series of PDFs, and including the PDF uncertainty, we can estimate the 
uncertainties of cross section calculations. 

For any calculated quantity a, we denote the central prediction, corresponding to 
a s (Mz) = 0.118, by a . There are two contributions to the uncertainty. The symmetric 
PDF uncertainty, denoted Aopdf, is calculated from the Hessian error PDFs by the "master 
formula" [1] 



. 1 

AlTpDF — ~> 



Here d is the dimension of the parameter space, and af 1 is calculated with the eigenvector 
PDFs. 2 The a s uncertainty of a is 



Aa as = -\/[a (A^)-a (A 2 )] 2 . (19) 



The combined uncertainty Aa for CTEQ6.6+CTEQ6.6AS is 



(Aa) 2 = (Aa PDF ) 2 + (Aa a f . (20) 



For the CTEQ6.6FAS series, the full uncertainty is computed according to Eq. ([18]) for 
d = 23. Then the prediction for the quantity a is 

a = a ± Aa. (21) 

Table [T] lists the predictions for a sample set of NLO cross sections at the Tevatron 
and LHC energies. Four processes are calculated: ti production [l8|, 19[, using a program 



from Ref. [19j; and Standard Model Higgs boson production, gg — > H, for Higgs masses 



Mn — 120, 160 and 250 GeV [20|. The production cross sections are predicted for three LHC 
energies: 7, 10 and 14 TeV. The central predictions a§ and uncertainty ranges Aa are given 
first for CTEQ6.6+AS (fitting method 1) and then for CTEQ6.6FAS (fitting method 2). 
Clearly, the predictions by the two methods are very close, even though not identical due to 
secondary effects like deviations from the quadratic approximation. 

In fitting method 1, we can also compare the relative sizes of the PDF uncertainty and the 
a s uncertainty. For processes dominated by the gluon scattering, such as the ti production 
or Higgs production examined here, the a s and PDF uncertainties can be comparable, as is 
observed in the Table. These results are illustrated for ti production at the LHC with energies 
7 TeV and 14 TeV in Fig. El and for Higgs boson production (Mh = 120 GeV) at the LHC 
with energies 7 TeV and 14 TeV in Fig. [71 The figures show the cross sections from individual 
CTEQ6.6(AS) eigenvector sets, as well as the resulting PDF and a s uncertainties, versus 



2 For these PDFs, based on CTEQ6.6, the number of fitting parameters is d = 22. PDF sets +i and — i are 
variations of the central fit corresponding to displacements in the + and — directions along eigenvector i. 
The symmetric error formula, Eq. (|18[) . is sufficient for the comparison here; at other times, we would use 
more detailed asymmetric error estimates, Eq. (|17[) . 
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do ± Act 


Tevatron, 1.96 TeV 


0.055 


0.0099 


0.0044 


0.011 


0.058 ± 0.012 


LHC, 7 TeV 


2.30 


0.085 


0.081 


0.12 


2.32 ±0.12 


LHC, 10 TeV 


5.08 


0.14 


0.15 


0.21 


5.10 ±0.22 


LHC, 14 TeV 


10.03 


0.26 


0.27 


0.37 


10.04 ±0.41 



TABLE I: Predicted cross sections for it production and Higgs boson production at the Tevatron 
and LHC, for different scattering energies. Cross sections are given in picobarns. The renormal- 
ization and factorization scales are chosen to be the top quark mass and Higgs boson mass for the 
it and gg —> H cross sections, respectively. The columns present the central prediction o"o, PDF 
uncertainty Apup, a s uncertainty A Qs , and total PDF±a s uncertainty Aa, computed according 
to the CTEQ6.6AS and CTEQ6.6FAS methods. 

the corresponding a s {Mz) values. The overall prediction based on the full CTEQ6.6+AS 
eigenvector set, <r ± Aa, is shown by an error bar: the inner bar is the PDF error alone; 
the outer bar is the combined PDF + a s error. According to the figures, the a s uncertainty 
of the total cross section in the shown processes constitutes between 70% and 110% of the 
PDF uncertainty. 
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FIG. 6: Cross section for tt production at the LHC, at center-of-mass energies of 7 TeV and 14 
TeV, as a function of a s (Mz)- The predictions for the the central set and the 44 eigenvector PDF 
sets are shown for a s {M z ) = 0.118. For the other values of a s {M z ) (0.116, 0.117, 0.119, 0.120), 
only the central prediction is shown. The combined uncertainty range (CTEQ6.6+CTEQ6.6AS) is 
shown as the error bar; cf. Table HI The inner error bar is the PDF error alone, and the outer error 
bar is the combined PDF + a s uncertainty. 

IV. CORRELATION BETWEEN a s AND THE PDFS 

The independence of the a s uncertainty from the PDF uncertainty in the CTEQ6.6AS 
method does not preclude existence of some correlation between the a s and PDF parameters. 
This correlation arises from the hadronic scattering experiments, which probe a variety 
of combinations of the PDFs and a s . Let us now examine which PDF flavors are most 
affected by variations in a B (Mz), and which scattering experiments impose the most relevant 
constraints. 

For this purpose, we go back to the full fit CTEQ6.6FAS with the floating a s and examine 
the correlation between a s (Mz) and individual PDF f a (x,Q) using the method outlined in 



17 




10.2 



0.116 



0.117 



0.118 



a, 



0.119 



0.12 




FIG. 7: Same as Fig.[6l for production of Standard Model Higgs boson with mass 120 GeV at the 
LHC. 



Refs. pi QUI- Let y} ±] (for % = 1, 23) denote one of 46 eigenvector PDFs f a (x, Q) for a 
chosen a, x, and Q, corresponding to the maximal acceptable displacements of orthonormal 
PDF parameters {zi, ...,223} from their best-fit values in the positive (+) and negative (— ) 
directions, respectively. Given 46 values of a s (Mz) (denoted as X^), we compute the 
correlation cosine, 



23 



cosy? 



4AX AY 



i=i 



•(+) 



X, 



(-) 



Y; 



( + ) 



Y; 



(-) 



(22) 



where AX and AY are the symmetric PDF errors, found from Eq. (1181) for d = 23. Values of 
cos ip close to +1, 0, —1 indicate strong correlation, no correlation, and strong anti-correlation 
between a s (M z ) and the PDF f a (x,Q) in question. 

Fig. [8] shows cosip versus x, for the PDFs that have the largest correlations with a s (Mz), 
at Q = 2 and 85 GeV. Other PDFs have little correlation with a s (Mz) and are not shown 
in the figure. We also identify the experimental data sets that cause the largest correlations. 
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CTEQ6.6FAS: correlation o\a s (M z ) with f~ a (x,Q) 
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FIG. 8: Correlation cosine cose/? between the best-fit QCD coupling a s (Mz) and PDFs f a (x,Q), 
plotted as a function of the momentum fraction x at Q = 2 GeV and Q = 85 GeV. 



As shown in Fig. [HI at Q — 2 GeV , the most significant correlation (or anticorrelation) of 
a s (M z ) occurs with: 

• the gluon PDF g(x,Q) (anticorrelated at x ~ 0.01 due to neutral-current DIS con- 
straints from HERA); 

• the singlet PDF T,(x, Q) (correlated at x ~ 0.4 due to constraints imposed by BCDMS 
and NMC neutral-current deep- inelastic scattering (DIS) data); 

• the heavy-quark PDFs, c(x, Q) and b(x, Q) (correlated at x = 0.05 — 0.2 due to con- 
straints by HERA charm and bottom semi-inclusive DIS data). 

At Q = 85 GeV, these correlations are reduced by the PDF evolution, with the exception of 
the heavy-quark PDFs. The same conclusions can be drawn from the simpler CTEQ6.6AS 



analysis, as well as from the correlation analysis by NNPDF [ll| . 

The best-fit value of a s (M z ) is thus determined by several types of the data, probing the 
gluon evolution in DIS at moderately small x, the singlet PDF evolution in DIS at large x, 
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and HERA charm semi-inclusive DIS data. The correlation of each kind disappears if the 
relevant data set is removed. For example, it is believed that the low-Q/large-x BCDMS and 
NMC DIS data prefer somewhat lower a s (Mz) values than the rest of the experiments (2l| . 
indicating the possible presence of higher-twist terms that are not explicitly included in most 
PDF analyses (22|. With the most suspect part of these data excluded, the central a s (Mz) 
value indeed increases to 0.119 - 0.120; the spike in cosip for the singlet PDF S(x, Q) at 
x ~ 0.4 also disappears from Fig. [HI 



V. CONCLUSION 

We conclude with this important point: while there are correlations between a s (Mz) and 
some PDFs, as demonstrated in Figure [HJ the total uncertainty of the CTEQ6.6FAS fit is 
essentially the same as the CTEQ6.6+CTEQ6.6AS fit, as demonstrated in Figure The 
CTEQ6.6+CTEQ6.6AS method captures all the correlation of a s and PDF parameters as a 
mathematical consequence of the quadratic approximation, upon which the Hessian method 
is based. The ability of the quadrature method to reproduce the total PDF+a s uncertainty 
(also observed in a related NNPDF study |Tl[) is thus more than a coincidence. 

We may predict cross sections using the CTEQ6.6FAS PDFs, with their 23 x 2 Hessian 
eigenvector PDF sets. The results are shown in the final column of Table [H The predictions 
agree well with those of the simpler method based on the CTEQ6.6+CTEQ6.6AS PDF 
sets. This again justifies our recommendation to use the CTEQ6.6+CTEQ6.6AS PDFs to 
estimate uncertainties induced by the PDFs and a s , according to Eq. fl2TJ|) . 

The experimental measurements included in the world average give a smaller a s uncer- 
tainty than the hadronic scattering data can match. Therefore, when making predictions for 
the LHC or Tevatron colliders, the most realistic estimate of the a s uncertainty will result 
from the CTEQ6.6+CTEQ6.6AS a s series, which covers 

a s {M z ) = 0.116, 0.117, 0.118, 0.119, and 0.120. (23) 

The central value of this series (assumed in the CTEQ6.6 eigenvector set) is closest to the 
world average. The PDFs for a s values at 0.116 and 0.120 provide an estimate of the a s 
uncertainty at approximately 90% CL, which is to be combined in quadrature with the 
CTEQ6.6 PDF uncertainty. The CTEQ6.6AS PDF sets are available from the CTEQ6.6 
PDF website ^ and as a part of the LHAPDF library Q. 
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Appendix A: Error matrix diagonalization for the QCD coupling 



In this Appendix we provide a proof of the method advocated in this study, whereby an 
uncertainty induced by the small variation of a s value is added in quadrature with the fixed- 
a s PDF uncertainty. The only necessary assumptions are the usual quadratic approximation 
for x 2 , and the understanding that the PDF sets for the "up" and "down" variations of a s 
are obtained through fits in which all of the PDF parameters are free. 

Assume we have an initial fit in which there are N free parameters, with an additional 
parameter held fixed. The additional parameter can be any additional degree of freedom 
that one wishes to study. Here it is intended to be a s [Mz)\ but no special features of that 
parameter will be invoked. 

We begin by carrying out the traditional Hessian diagonalization procedure, applicable 
in the vicinity of the best-fit (BF) combination of the PDF parameters. Starting from the 
"shape" parameters that control the input PDFs at scale Qq, we define new fitting parameters 
ai, . . . , ajv that are coefficients of the eigenvectors of the original Hessian matrix, to obtain 



X 



Xbf 



N 

i=i 



(Al) 



We now wish to include an additional degree of freedom in the fit. Hence we introduce a 
new parameter do, which could be defined as proportional to the deviation of a s (Mz) from 
its best-fit value, a s (Mz) — a s> BF(^z); or better, as proportional to hx[a s (Mz) / ch 8 ,bf{Mz)]- 
Including the new degree of freedom, we have, in general, 



N N 

X 2 = Xbf + ba o + ^^Hijaidj , 

i=0 j=0 



(A2) 



where 



H 



Here 



/l 

Pi 1 

P2 



Pi P2 P3 



1 



P3 1 



\Pn o 



Pi 



Pn\ 






1/ 



(A3) 



(A4) 



H describes the correlation between ao and aj, as it gives the variation of \ 2 i n t ne { a o, a i}~ 
subspace. The linear term bao in flA2j) allows for the possibility that the new minimum 
might not be at a = 0. Note that b = (dx 2 /da ) BF , which describes how x 2 varies with a 
at the original best fit. The new diagonal element H Q0 was chosen to be 1 by including an 
appropriate scaling factor in the definition of a . 
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As do varies, the PDF parameters a±, are adjusted so that to reduce the increase in 
X 2 ■ The direction along which this increase is minimal (direction AC in Fig.[TJ corresponds 
to dx 2 /dcii = for % — 1, ...N, i.e., it is a line given by = —pia according to Eqs. (1A2I) 
and flMD. 

Define new variables 

Zi = ai+^ao, for i = 1, N, (A5) 

a + ; b (A6) 
2 VI -C 2 



where 



C 



N 

\ ( A7 ) 

Since \ 2 is assumed to have a minimum, hence, < C < 1. The transformation (1A51 IA6[) to 
the Zj coordinates consists of an ao-dependent translation (shear) along the PDF coordinates 
aj, followed by a scaling and a translation along the a s coordinate clq. In the new coordinates, 
X 2 is diagonal, 

b 2 N 

X 2 = Xbf " m _ C 2) + E ^ ( A8 ) 

^ ' i=0 

Variations along zq and ^ for i = 1,...N are explicitly independent, and we can calculate 
the symmetric uncertainty AX for an observable X(zq, Zi, zn) by 

1 N 

(AX) 2 = - [^(0, Zi = T, 0) - X(0, ..., Zi = —T, 0)] 2 , (A9) 

j=0 

for tolerance Ax 2 = T 2 . 

Below, we will focus on the case that 6 = 0, which is valid in this paper. From Eqs. ( 1A5I) 
and ( 1A6|) . the extreme values of ai consistent with A% 2 = T 2 are af = zf = ±T. They occur 
at a = zq = (i.e., at the best-fit a s ) and Zk = for ^ i. The extreme values of a occur 
at ^ = ±T and ^ = for i = 1, ...,N, i.e., Oq = ±T/\/T — C 2 and a^ 1 = — PiO^ . Hence 
Eq. (1A9P corresponds exactly to adding in quadrature the uncertainty based on the allowed 
range of computed at a = (analogous to the CTEQ6.6 uncertainty in our application) 
to the uncertainty based on the allowed range of ao, computed with at their preferred 
values for the extreme excursions in a (analogous to the CTEQ6.6AS uncertainty). 

An alternative way to construct the Zi basis is to find the eigenvectors of the (X+l) x (N + 
1) Hessian matrix H in Eq. (1A3I) and then perform an additional orthogonal transformation. 
To find the eigenvectors, note that H can be written in a block form 

h -(p P i T ) < A10 > 

where 1 is a unit element, p T is a row vector (pi,p 2 , ...,Pn), and I is an N x X unit matrix. 



One eigenvector is 

V (°) = F ( 

P 

corresponding to an eigenvalue Aq = 1 + C. The normalization factor is F = l/(v / 2C). 



Vi°) = F ( C ), (All) 
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V (l) = F l~ c \ (A12 ) 



A second eigenvector is 

V« — F ( 

P 

with an eigenvalue Ai = 1 — C. 
• Now consider vectors of the form 

V=(°) (A!3) 

where p T q = 0. These are orthogonal to the first two eigenvectors; the iV-dimensional 
vector q is orthogonal to p. For any q, we have H V = V. All their eigenvalues are 
A J = l,forj = 2,3,...,iV. 

So we have a set of eigenvectors of H, denoted by V^°\ V^, 'V <y2 \ . . . , V^, and satisfying 
the orthonormality condition 

N 



fc=0 



The next step is to use the eigenvectors to define new coordinates to, . . . , tjy by 

N 

a* = E *i v ^ ■ ( A15 ) 

3=0 

It is easy to compute \ 2 from Eq. (IA20 by making use of the eigenvector properties 

HV (J) = Xj V {j) , (A16) 

where Ao = 1 + C, Ai = 1 — C, and A^ = 1 for j = 2, . . . , N; and then using the orthogonality 
property (lA14j) . The result has the expected diagonal form 

N 

X 2 = xl F + (1 + C% 2 + (1 - C)tl + J^ti > (A17) 
which is further simplified to 

N 

X 2 = X 2 BF + ( A18 ) 
i=0 

by rescaling of the first two coordinates: 

= Vl + Ct ; y 1 = Vl-Ct 1 ; Vi = U fori = 2,...,N . (A19) 

The function \ 2 is explicitly diagonal in the new coordinates yi, but not in the parameter 
cto that is of interest here, given that 



a = 2/ / V2(l + C) - yi /V2(l - C) . (A20) 
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It can be made so by a further orthogonal transformation called the Data Set Diago- 
nalization [l7[. The general method for finding the necessary transformation is to write 
cXq = J2i=o X^=o Wi Wj and use the eigenvectors of G as new basis vectors. However, in 
the present case, the required transformation can be found more simply by observing that 
one of the new coordinates must be proportional to do, and the others must be orthogonal 
to it. The result is 

z o = VWvo - VWvi, 
*i = VWvo + VWvii 

Zj = Vj for j = 2,...,N, (A21) 

which yields 

N 

X 2 = XI F + J2 Z 1 ( A22 ) 
3=0 

equivalent to Eq. ( ]A8j) with 6 = 0. Because \ 2 is a sum °f independent terms for each Zi, 
the uncertainties of these variables should be combined in quadrature. This conclusion does 
not depend on the value of C, i.e., the strength of correlation between the new degree of 
freedom ao and the original parameter space of Oj in Eq. flAll) . The dependence of x 2 on 
ao comes entirely from the coordinate Zq. The space spanned by Z\, . . . , at Zq — has 
a = 0, so it is the same fitting space spanned by ai, . . . , in the original iV-parameter fit. 
Hence the uncertainty associated with z±, . . . ,Zn is the same as the uncertainty associated 
with ai, . . . , Oat. It is therefore not necessary to explicitly carry out the transformations, and 
the overall formula for the uncertainty is the result we have used in Eq. (I2UD . 
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